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Abstract- In this paper, two robust controllers are designed for the 
synchronization of Lorenz-Rossler chaotic systems. First an active 
sliding mode control is presented in which the resulted control 
parameters not only secure the stability of the error system; they 
can also be adjusted to obtain the desired rate of error 
convergence. Second, linear quadratic gaussian/loop transfer 
recovery (LQG/LTR) controller is proposed which allows the 
designer to shape the principal gains of the return ratio at either 
the input or the output of the plant, to achieve robust stability and 
robust performance specifications. The numerical simulation 
results show the effectiveness of two proposed methods. 
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I. INTRODUCTION 

Chaos is an interesting complex dynamical phenomenon 
which plays an important role in the field of non-linear science. 
Recently, the traditional trend of analyzing and understanding 
chaos has evolved into a new phase. Research in this field has 
been moved to chaos control, synchronization, modeling, and 
even chaotification [1-7]. 

Synchronizing the chaotic systems has gained a lot of 
attention among scientists from a variety of research fields in 
recent years [1-7]. The chaos synchronization is essential in 
some applications such as communication and secure 
communication [4]. The idea of synchronizing two identical 
chaotic systems was first introduced by Carroll and Pecora [1, 
8]. A wide variety of approaches such as active control [9-13] 
and sliding mode control [15, 16] have been employed for the 
synchronization of chaotic systems, most of which are 
designed to synchronize two identical chaotic systems. Bai and 
Lonngren used active control technique to synchronize two 
Lorenz systems [9]. The same approach, because of its inherent 
advantages, has been used by many other researchers to 
synchronize other chaotic systems [10-13]. The sliding mode is 
another efficient method for synchronization [15, 16]. This 
method is a discontinuous control strategy that involves 
selecting a switching surface for the desired dynamics and then 
designing a discontinuous control law. The system trajectory 
first reaches the surface and then stays on it as long as there are 
no perturbations [14]. 

In this paper, the goal is to synchronize two different 
chaotic systems: Lorenz system as master and Rossler system 
as slave. To achieve the goal, two robust controllers are 



designed: Active sliding mode and LQG/LTR controllers. First, 
the active sliding mode controller is designed as a combination 
of two design procedures: the active controller and the sliding 
mode controller [17]. The technique requires two stages. The 
first stage is to select an appropriate active controller to 
facilitate the design of the sequent sliding mode controller. The 
second stage is to design a sliding mode controller to achieve 
the synchronization. In this paper, the active sliding mode is 
implemented and in order to improve the method introduced in 
[17] uncertainty and disturbance are added to the slave system. 
Secondly, LQG/LTR controller is designed to achieve the 
required specifications which are described as good tracking 
and disturbance rejection by augmenting the open-loop system 
with integral action in each loop, robust stability, and robust 
performance of the closed-loop system with respect to the 
uncertainty of the plant. The LQG/LTR method is rooted in 
optimal control theory and is a systematic design approach 
based on shaping and recovering open-loop singular values 
[18]. LQG design approach is the combination of the optimal 
state-feedback regulator (LQR) and the Kalman filter (KF). 
Then, to recover the properties of both methods, the loop 
transfer recovery (LTR) is utilized. The LQG approach 
accompanied by LTR technique is named LQG/LTR method 
[18]. 

The organization of this paper is as follows: Section II 
focuses on the description of the dynamics of the Lorenz and 
Rossler chaotic systems, while the introduction of the active 
sliding mode method and the simulation resultants of the 
synchronization of two chaotic systems are provided in Section 
III. The design of the LQG/LTR method and the corresponding 
simulation results are described in Section IV. The comparison 
of the simulation results of two methods is given in Section V. 
Finally, conclusions are presented in Section VI. 

II. THE DYNAMICS OF LORENZ AND ROSSLER 
CHAOTIC SYSTEMS 

The Lorenz system is known as a simplified model of 
several physical systems. The Lorenz equations may describe 
such different systems as laser devices, disk dynamos and 
several other problems related to the convection [19]. The 
system is the first chaotic system which was introduced by 
Edward Lorenz in 1963. The Lorenz system is described by the 
following nonlinear differential equations: 
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(1) 



x = a(y-x) 

y = -xz + yx-y 
z-xy-fiz 

where ' " and z are state variables and ' " and ' are three 
real constants. Figure 1 illustrates that the system has chaotic 



attractor when 



a=10,/? = 8/3 and 7 = 28 
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(d) 

Fig. 1 Lorenz chaotic attractor: a) projection on x-y, b) projection on y-z, c) 

projection on x-z, d) in x-y-z space 

In the 1970s, Otto Rossler, a non-practicing medical doctor, 

became interested in chaos and concocted a three-dimensional 

autonomous system even simpler than Lorenz [20]. The 

dynamics of the Rossler system is described by: 



x=-y-z 
y = -x + ay 
z = b + z(x - c) 



(2) 



where ' ^ and z are state variables and ' and c are system 
parameters. Figure 2 shows the trajectories of the Rossler 
attractor where a = o = 0.2 anc j c = 5.7 

III. ACTIVE SLIDING MODE CONTROLLER DESIGN 

The active sliding mode controller is a combination of two 
design procedures: the active controller and the sliding mode 
controller, which are given in succession. 

A. Active Sliding Mode Controller Design 

Consider the master chaotic system described by the 
following nonlinear differential equation: 

x(t)=f 1 (x) (3) 

where *■ ' denotes the system's three-dimensional state 

vector and (1 ' is the system's nonlinear equation. 

The slave system is described as below, in which *■ ' as 

the control input and *■ ' *' e as the uncertainty and 
disturbance term are added to it: 



y(t) = f 2 (y)+u(t) + M(t,y) 

>3 



(4) 



where " ^ ' is the slave system's three-dimensional state 

f ' R^ —> R^ f 

vector and ' 2 ' imply the similar role as '* for the 

master system. If ' l ~ ' 2 , then x and " are the states of two 
identical chaotic systems. Otherwise, they represent the states 
of two different chaotic systems. We consider that: 



M(t,y) = Af(y) + d(t) 



(5) 



where Af e J? anc j •. > are the uncertainty and the 

disturbance of the slave system, respectively. It is considered 
that the uncertainty and disturbance are bounded, such that: 



M <M i 

\\d(t)\\<M 2 



(6) 

(7) 



where 1 ' 2 are the vectors of constant scalars. Then 

we can assume that: 



0<M 1 +M 2 <M 3 



(8) 



where 3 is a constant vector. 




(a) 
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(b) 




(c) 




-20 -10 



(d) 



Fig. 2 Rossler chaotic attractor , a) projection on x-y, b) projection on y-z, c) 
projection on x-z, d) in x-y-z space 



Therefore, we have: 



M (t,y)\\<M 3 . 



(9) 



As y'J ' represents the disturbance and uncertainty, 
which is unknown, but according to (9), it is bounded and its 



error dynamics [17]. In other words, the input vector is 
considered as: 

u(t) = Ae(t) + £(x) - f 2 (y) + Kw(t) (12) 

where A is a canonical matrix and represents the linear part of 

>u a ■ K =[k u k 7 ,kof sR 3 . t t . 

the error dynamics, <- l, i, 51 ls a cons t an t g ain 

vector and *■ ' is the control input that satisfies: 



W(t): 



w + (t) S(e)>0 
w"(t) S(e)<0 



(13) 



and ^ ' is a switching surface which determines the 

desired dynamics. 

The error system (10) is then rewritten as below: 

e = Ae(t) + Kw(t) + M(t,y). (14) 

In next section, the appropriate sliding mode controller will 
be designed according to the sliding mode control theory. 

1) Sliding Surface Design 
The sliding surface can be defined as follows: 



S(e) = Ce 



(15) 



C T =[c 1 ,c 2 ,c- i f sR 3 



where l' 2 ' 3 is a constant vector. The 

equivalent control approach is found by the fact that *• ' ~~ 
is a necessary condition for the state trajectory to stay on the 

switching surface ^ ' ~ . Hence, in the sliding mode, the 
controlled system satisfies the following conditions: 



and 



S(e) = 

S(e) = 0. 
Now, using (14), (15) and (17), one can obtain 



(16) 
(17) 



S(e) = 



dS(e) 



, e\t)=C[Ae(t) + K W (t) + M(t,y)] 

oe (18) 



= 0- 

Solving (18) for ^ ' results in the equivalent control 

w eq(t) , , 
H , as below: 



norm is less than 



M- 



The synchronization problem is to design a controller 
which synchronizes the states of the master and slave systems. 
The dynamics of the synchronization error can be expressed as: 



w(t) = -(CK)\CAe). 



(19) 



e=f 2 (y)-f 1 (x) + u(t) + M(t,y) 



Replacing eq for ^ ' in (14), the state equation in 
the sliding mode is as follows: 



(10) 



where 



■ y -x 



e(t) = [I-K(CK) 1 C]Ae 



(20) 



The goal is to design the controller *■ ' such that: 

Iim|e(t)| = . (li) 

According to the active control design procedure, one uses 



the control input 



U(t) 



to eliminate the nonlinear part of the 



As long as the system (20) has all eigenvalues with 
negative real parts, it is asymptotically stable. For this purpose, 
the matrices K and ^ must be chosen properly. 

2) Design of the Sliding Mode Controller 

We assume that the constant plus proportional rate plus 
uncertainty norm reaching law is applied. The reaching law can 
be chosen as below: 
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S(t) = -CK sgn(S ) - rS +CM 3 
where ° *•'' 



(21) 



M 3 eR" 



denotes the sign function and 3 with 

positive constants is the additive norm of the uncertainty and 

disturbance term. The gain vectors ' and the scalar r 
must be chosen such that the sliding condition is satisfied and 
the sliding mode motion occurred. 



From (14) and (15), it can be found that 
S (t) = C[Ae(t) + Kw (t) + M (t ,y)]. 



(22) 



Now, from (21) and (22), the control input is determined as 
below: 

w (0 = -{CK)\C (rl + A)e(t) + CK sgn(S)]- (23) 



B. Stability Analyse 

The stability of the proposed method is studied in two parts: 
the convergence to the switching surface and the stability of the 
closed-loop system, which are given in the following. 

1) Analyze of the Convergence to the Switching Surface 

The error system in the sliding mode is determined as 
follows: 



e(t) = [A-K(CKy 1 C(rI + A)]e(t) 
-Ksgn(S) + M(t,y) 



(24) 



which is a three-dimensional linear system. The error dynamics 
can be reduced to a two-dimensional one, as: 



-I< 



i=l 



(25) 



C • g — y — v . 

where ! is a positive scalar and ' J ' ! 
We suppose that 

|-Ksgn(S)+M(t,y)|<0 



(26) 



,3 . 



where $eR is a constant vector with positive scalars. 



Substituting (25) and (26) into the system (24), one can 
obtain: 



E = BE + O 



(27) 



where £ =[e 1 ,e 2 ] e R and B e R x . Using sub-matrix 
multiplication, the i,j th entry b ; . of matrix B can be 
calculated easily: 

bij =Oi,i -o i>n Cj -{CKy'K^Cai -c } (Ca")] 



i,J=l,2, 



n 



where IJ is the i,j th entry of the matrix ■"■ and a denotes 

the j th column in -A . All the eigenvalues of the system (27) 
has negative real parts. So, all the states of the system converge 
to the sliding surface. 

2) Analyze of the Stability of the System 



To check the stability of the controlled system, we consider 
the following Lyapunov candidate function: 



v = -s 2 . 

2 
The time derivative of (29) is as follows: 

V =SS = S (C (Ae(t) + Kw (t) + M (t,y))) . 
By replacing (23) into (30), we obtain: 



V = -CKS sgn(S )-rS 
By choosing 

CK>CM 3 -rS 



-SCM- 



(29) 



(30) 



(31) 



we have v < (J The inequality (32) can be simplified to: 



CK >CM. 



(32) 



(33) 



The matrices ^ and K must be chosen such that (33) 
holds. According to the Lyapunov theory, the above property 

implies the boundedness of the sliding surface ^ . 

C. Simulation Results 

Active sliding mode approach is applied to synchronize 
Lorenz and Rossler chaotic systems. The technique capability 
is examined by numerical simulations, which is discussed in 
the following. 

To observe the synchronization behavior in the Lorenz and 
Rossler systems, we assume that the Lorenz system drives the 
Rossler system which contains uncertainty and disturbance. 
Therefore, the master and slave systems are defined, 
respectively, as follows: 



x 1 - a(x 2 -x x ) 

2 = -x ]X 3 + yx 1 - X 2 



(34) 



X 1 X 2 



fix. 



and 



yi = -y 2 -y3+ d i(0+Af 1 (y)+u 1 (t) 
y2=yi+a> , 2+ d 2(0+Af 2 (y)+u 2 (t) 
ys^+ysOi-O+^co+AfgOo+i^t)- 



(35) 



By subtracting (34) from (35), we obtain the error dynamics 
as below: 



?i = -y 2 -.y3+a(x2-Xi) 

+ d 1 (t) + Af 1 (y) + u 1 (t) 



(36) 



+ d 2 (t) + Af 2 (y) + u 2 (t) 

e 3 =b + y 3 (y 1 - c) + x x x 2 + j3x 3 

+ d 3 (t) + Af 3 (y) + u 3 (t) 

where d i (t) is external disturbance which is assumed to be a 
Gaussian white noise and Af ; is uncertainty which is 
considered as below: 
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Af 1 (y) = 0.5jjsin(,>y,.) 
i=l 

3 

Af 2 (y) = 0.5jJsln(J^ / ) 

3 

Af 3 (y )=0.5jjsln(f^,) 



(37) 



"o 


1 


o" 








1 












(38) 



The control parameters are chosen as below: 

0" 

1 



K = [0 1 if 

C =[0.5 0.5 0.1] 

r =20- 

As it is mentioned before, the matrices K and ^ are 
chosen such that the inequality (33) holds. 

The above selections result in negative eigenvalues for the 

error dynamics. The control signal *• ' is determined as 
follows: 

w(0 = [-16.6667 -17.5 -4.16667]e(t) + sgn(S (e))- (39) 

We selected the parameters of the Lorenz system as 

a=10,B = 8/3,y = 28 A ., . f ., D .. , 

' r " and the parameters of the Rossler 

system as a = t> = U.z anc j c = b./ _ j ne initial values of the 

master and slave systems are 

X 1 (0) = 0,X 2 (0) = -0.M,X 3 (0) = 9 and y 1 (0) = -9,y 2 (0) = > 

J3y)~ t respectively. These choices result in the initial error 
valuesof e 1 (0) = -9 ) e 2 (0) = 0.01, e 3(0) = -9_ 




10 15 20 

time(sec) 



Fig. 3 Control signal w (t) where controller is activated at t=l sec 




time(sec) 

Fig. 4 Active sliding mode switching surface where controller is activated at 

t=l sec 




Fig. 5 Synchronization error signals where active sliding mode controller is 
activated at t=l sec 
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10 20 30 40 50 




20 30 

time (sec) 



Fig. 6 Active sliding mode control signals applied at r=l sec 

The simulation results are given above. Figure 3 shows the 
control signal w (t) which is applied at t=l sec, while the 
switching surface is depicted in Figure 4. The error dynamics 
of the synchronization is illustrated in Figure 5. The control 
signals u (t), i =1,2,3 which are applied to the Rossler 

system are shown in Figure 6. 

IV. LQG/LTR DESIGN METHOD 

A. Problem Definition 

Suppose that we have a plant model in state space form, as 
follows: 

f X = Ax + Bu + Tw 
[y = Cx + Du + v 
where W and V are 'white noises', and 



(40) 
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E{ww T }=W >0 
E{vv T }=V >0 
E {wv } = • 



(41) 



The problem is to design a feedback control law, which 
minimizes the following cost function: 



J = lim£{ f (z T Qz + u T Ru)dt} 



(42) 



where z "^ x is some linear combination of the states and 

Q =Q r >0, R = R T >0 ., .,.. , . riQ1 

are the weighting matrices L18J. 

The solution to the LQG problem is presented by the 
separation principle, which states that the optimal result is 
achieved by the following procedure: first, one should obtain 

an optimal estimate x of the state x and then use this estimate 
as if it were an exact measurement of the state to solve a 
deterministic linear quadratic control problem. This procedure 
reduces the problem in two sub-problems, the solution to which 
is known [18]. 

B. The Solution of the LQG Problem 

As mentioned above, the separation principle reduces the 
problem in two sub-problems. The first sub-problem is the state 
estimation. It is addressed by Kalman filter theory. The second 
sub-problem is to find the control signal which will minimize 
the following cost function: 



\(z T Qz + u T Ru)dt. 



On the assumption that 



(43) 



(44) 



the solution to this is to let the control signal u be the linear 
function of the states, as below: 



u = -K c x 

K 



(45) 



where c is known as the state-feedback gain matrix [18, 21]. 
The optimal state-feedback matrix c is given by: 

K c = R 1 B T P c (46) 

where Pc satisfies the following algebraic Ricatti equation: 

A T P C +P C A-P C BR 1 B T P C +M T QM =0 (47) 

P = P T > K 

and c c . The Kalman filter gain matrix f is given, 

as below: 



K f =P f C T V~ 1 



(48) 



where ' satisfies the following algebraic Ricatti equation: 
P f A T +AP f -P f C T V 'CV x P f +TWT T = (49) 
P f =P/>0 r 

and f f [18]. 



The resulting closed loop system is robust against any 

unstructured multiplicative uncertainty A j at the input of the 
plant, for which we have: 



o-(A(;w)) 



< 



1 



(50) 



C. Loop Transfer Recovery (LTR) 

LQG design approach is the combination of the optimal 
state-feedback regulator (LQR) and the Kalman filter (KF) as 
illustrated in Figure 7. Since both the LQR and KF methods 
have such good properties, it might be expected that LQG 
compensator would generally yield good robustness and 
performance. Unfortunately this is not the case: as Kalman 
filter loop gives good robustness and tracking properties, so we 
are in search of a compensator that recovers all these properties 
at the plant input. Similarly optimal state-feedback controller 
exhibits good robustness in the presence of unstructured 
uncertainties, therefore it is necessary for a Kalman filter to 
recover all these properties at the plant output. This procedure 
is called Loop Transfer Recovery. The former is called loop 
transfer recovery at the plant input and the later is called LTR 
at the plant output. At this stage, the dynamics of the plant are 
inversed and that is why LQG/LTR procedure is not valid for 
non square plant. During this, some of the filter eigenvalues are 
placed at the zeros of the plant, the remainder being allowed to 
be arbitrary fast. 



Compensator 

K(s) 



r=0 _l + 



+ Al iL 



-m^®+ % *x»> 



i' 



++ i 



c 4 



Cf.sj 



Plant 



Fig. 7 Complete diagram of LQG/LTR method [9] 

D. General Design Procedure 

The following two step design approach is used for the 
LQG/LTR controller design for square plants. 

Step-1: Design a Kalman filter by manipulating the 

covariance matrices W and * until the return ratio 

-C(sZ-A) _1 2C f ,_ 

' obtained is satisfactory at the plant output. 

Emphasis should be put on the aspects such as 0-dB cross over 
frequencies for the principal gains, possibly balancing the 
principal gains [18]. 

Step-2: Design an optimal state-feedback by 



setting 



M=C,Q = Q +qI 



and 



R = pl 



(alternatively 



choose ' ). Increase " (or reduce " ) until the 

return ratio of the compensated plant has converged 

«, • ■ , , -C{jwl-A) l K f ff . . , 

sufficiently closely to ' over a sufficiently 

large range of frequencies [18]. 
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It is often desired to design the return ratio at the plant 
output instead of the plant input. Otherwise, optimal state 
feedback matrix has to be designed in step-1 and then Kalman 
filter in step-2. 

E. Simulation Results 

Based on the procedure described in the previous sub- 
section, an LQG/LTR controller is designed for the 
synchronization of the Lorenz-RSssler chaotic systems. Loop 
transfer recovery is designed at the plant output, it means that, 

K f 

first the Kalman filter matrix ( ' ) is calculated by setting 

— ' * — 1 and by shaping the principal gains in four 
steps as discussed in [18], then the optimal state feedback 

matrix ( c ) is obtained by setting ' ^ and 

" . In the Kalman filter design procedure, we manipulate 
the singular values of the open-loop plant. 

As the open-loop system is augmented with integral 
action in each loop, the closed-loop system has good tracking 
and disturbance attenuation according to simulation results. 
Figure 8 compares the principal gains of the return ratio of the 
original error system and the system after shaping the principal 
gains in four steps. 
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Fig. 8 Principal gains of the open-loop return ratio of the error dynamics; RRfl: 
original error system, RRf4: system after shaping the principal gains in four 

steps 

The resulting observer gain matrix and optimal state 
feedback gain matrix are as below: 



K f 



K, 



12.0741 


-7.4608 


0.2312 


-7.4608 


12.9679 


0.2631 


0.2312 


0.2631 


4.4504 


108.2113 


-101.9002 


1.0998 


-98.0712 


111.6867 


-2.2043 


0.4950 


0.6117 


9.9645 


1000 0.5 


1 


0" 




0.5 1000 


0.5 


1 




0.5 


1000 


1 





The designed controller in state space form 

[4,B C ,C C ,D C ] . fll 

c c c c is as follows: 



A =A-B,K -K f C+K f DK 
B =K f 
C =K 

c c 

D =0 



(52) 



r A R C D 1 

where L ^' 1' v lJ is the augmented system model. The 

order of the designed controller is 6. 

In this paper, LQG/LTR control method is implemented to 
synchronize two different chaotic systems. The designed 
controller is tested in MATLAB simulink, in which a Gaussian 
white noise is applied at the plant input. 

The error signals, which are shown in Figure 9, confirm 
that Rbssler system is synchronized with Lorenz system. 
Figure 10 depicts corresponding control signals which are 
applied at t=l sec. 




20- 




10 15 20 

time(sec) 

Fig. 9 Synchronization error signals where LQG/LTR controller is activated at 

t=l sec 
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10 20 30 40 



O | ^^vrnw^yiJ^/Jty^JflW^ 



10 20 30 40 



" O L. -^.i^^^fj^A^fl,^^^^^^ 



20 30 

time (sec) 



Fig. 10 LQG/LTR control signals applied at r=l sec 
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V. THE COMPARISON OF THE RESULTS 

In this paper, the goal was to synchronize two different 
chaotic systems: Lorenz system as master and Rossler system 
as slave. To achieve this, two robust controllers were designed: 
active sliding mode and LQG/LTR controllers. First, the active 
sliding mode controller was designed as a combination of two 
design procedures: the active controller and the sliding mode 
controller. In this method an external disturbance which was a 
Gaussian white noise and uncertainty term considered as sine 
function of the state variables were applied to the slave system. 
The simulation results showed that two systems are 
successfully synchronized. Secondly, LQG/LTR controller was 
designed which was a systematic design approach based on 
shaping and recovering open-loop singular values and provided 
stability automatically. The corresponding simulations 
illustrated that the synchronization errors converged to zero 
very fast while a Gaussian white noise was applied at the plant 
input. 

As shown in Figure 11, the comparison of the error signals 
of two proposed methods demonstrates that the LQG/LTR 
control method provides better convergence rate to zero and 
less mean value of the error. According to Figure 11, the 

p p Q 

synchronization error signals 15 2 and 3 of the LQG/LTR 
control method reach zero when' — b.DSeC ; bsec anc j > bsec ( 

respectively, while the error signals v 2 and 3 of the active 
sliding mode method converge to zero when t—lZ.oSBC^ 
loSec anc j 11) sec ^ respectively. In addition, the comparison 
of the peak value of the error signals after applying control 

signals at ' — ^ sec confirms the advantage of the LQG/LTR 
technique. 

As shown in Figure 12, control efforts of both controllers 
for synchronizing two chaotic systems are in acceptable range, 
but the active sliding mode-based controller has chattering, as 
its inherent property, consequently LQG/LTR controller 
provides lower energy control signals. Besides, as both master 
and slave systems are chaotic, the control signals to 
synchronize the two systems have chaotic-like behaviour but 
they are bounded for all the times. 



VI. CONCLUSION 

In this paper, two robust controllers were used for the 
synchronization of Lorenz-Rbssler chaotic systems based on 
active sliding mode and LQG/LTR methods. The active control 
scheme had been successfully applied to synchronize both the 
response Rossler system and the drive Lorenz system, which 
provided efficient control law to make system robust with 
respect to system unmodeled dynamics and external 
disturbances. 
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Fig. 12 The synchronization control signals of two designed controllers; (solid: 
LQG/LTR method, dashed: Active sliding mode) 

The LQG/LTR controller was designed in two steps: first 
Kalman filter matrix was obtained by shaping the principal 
gains, and then the optimal state-feedback matrix was 
calculated. Due to the simulation results, the synchronization of 
Lorenz-Rbssler systems was achieved successfully using this 
method. As LQG/LTR method provided better convergence 
rate to zero with less mean value of the error signals, and lower 
energy control signals, this technique has more advantages. 
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Fig. 11 The synchronization error signals of two designed controllers; (solid: 
LQG/LTR method, dashed: Active sliding mode) 
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